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Abstract 



New nondiagonal G2 inhomogeneous cosmological solutions are presented 
in a wide range of scalar-tensor theories with a stiff perfect fluid as a matter 
source. The solutions have no big-bang singularity or any other curvature 
singularities. The dilaton field and the fluid energy density are regular every- 
where, too. The geodesic completeness of the solutions is investigated. 
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All versions of string theory and higher dimensional gravity theories predict the existence 
of the dilaton field which determines the gravitational "constant" as a variable quantity. 
The existence of a scalar partner of the tensor graviton may have a serious influence on the 
space-time structure and important consequences for cosmology and astrophysics. A large 
amount of research has been done in order to unveil the possible cosmological significance 
of the dilaton [|l|, 0- |jl^(and references therein). With a few exceptions most of the 
cosmological studies within the scalar-tensor theories were devoted to the homogeneous 
case. The homogeneous models are good approximations of the present universe. There is, 
however, no reason to assume that such a regular expansion is also suitable for a description 
of the early universe. Moreover, as is well known, the present universe is not exactly spacially 
homogeneous. That is why it is necessary to study inhomogeneous cosmological models. 
They allow us to investigate a number of long standing questions regarding the occurrence of 
singularities, the behaviour of the solutions in the vicinity of a singularity and the possibility 
of our universe arising from generic initial data. 

In this work we shall address the question of the occurrence of singularities in inhomo- 
geneous cosmologies within the framework of scalar-tensor theories. 

As is well known, most of the homogeneous models (both in general relativity and in 
scalar-tensor theories) predict a universal space-like big-bang singularity in a finite past. It 
was, therefore, believed that this would be the usual singularity in general. The inclusion 
of inhomogeneities drastically changes this point of view. There are inhomogeneous cos- 
mological solutions in general relativity which have no bing-bang or any other curvature 



singularity. The first such solution was discovered by Senovilla in 1990 Ref. Senovilla's 



solution represents a cylindrically symmetric universe filled with radiation. This solution 



has a diagonal metric and is also globally hyperbolic and geodesically complete |jT6[. Sen- 
ovilla's solution was generalized by Ruiz and Senovilla in Ref. [|T^ where a large family of 
singularity-free diagonal G2 inhomogeneous perfect fluid solutions was found. Nonsingular 
diagonal inhomegeneous solutions in general relativity describing cylindrically symmetric 



universes filled with stiff perfect fluid were found by Patel and Dadhich in Ref. [Iq]. Other 



examples of diagonal nonsingular solutions in general relativity can be found in Refs. 



In Ref. , Mars found the first nondiagonal G2 inhomogeneous cosmological solution of 
the Einstein equations with stiff perfect fluid as a source. This solution is globally hyperbolic 
and geodesically complete. Mars's solution was generalized by Griffiths and Bicak in Ref. 

Within the framework of scalar-tensor theories there are also inhomogeneous cosmological 



solutions without big-bang or any other curvature singularity. In Ref. PH|, Giovannini de 



rived gravi-dilaton inhomogeneous cosmological solutions with everywhere regular curvature 
invariants and bounded dilaton in tree-level dilaton driven models. In a subsequent paper 
p9| , it was shown that these solutions describe singularity-free dilaton driven cosmologies. 
A nondiagonal inhomogeneous cosmological solution with regular curvature invariants and 
unbounded dilaton in the tree level effective string models was found by Pimentel ||3D| . Very 
recently, inhomogeneous cosmological solutions without any curvature singularities were ob- 
tained by the author in a wide class of scalar-tensor theories with stiff perfect fluid as a 



source 



In this work we take a further step upwards and present new nondiagonal G2 inhomo- 
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geneous cosmological stiff perfect fluid solutions with no curvature singularities in a wide 
range of scalar-tensor theories. 

Scalar-tensor theories (without a cosmological potential) are described by the following 



action in Jordan (string) frame ||32|, 



S = I d'xy^ {F{^)R - Z{^)g'^''d,^d,^) (1) 

Here, is the bare gravitational constant and R is the Ricci scalar curvature with 
respect to the space-time metric 5^^,^. The dynamics of the scalar field $ depends on the 
functions -F($) and Z{^). In order for the gravitons to carry positive energy the func- 
tion -F($) must be positive. The nonnegativity of the energy of the dilaton requires that 
2F($)Z($) + 3[dF($)/d<l>]^ > 0. The action of matter depends on the material fields 
and the space-time metric g^^ but does not involve the scalar field $ in order for the weak 
equivalence principle to be satisfied. 

As a matter source we consider a stiff perfect fluid with equation of state p = p. 

The general form of the solutions is given by 

ds^ = F-\(^{t)) [e^""''" cosh{2at)i-dt^ + dr^) 

+ cosh(2at)rf0^ H -{dz + ar'^dcj)) 

cosh(2at) 



cosh(2at) 

= F-V2($(t))e(i/2)^'^'^' cosh'/\2at) 

The solution depends on three parameters - a, 7 (7 > 1) and A. The range of the 
coordinates is 

- 00 < t, 2; < 00, < r < 00, < < 27r. (3) 

The explicit form of the functions $(t) and /(A), and the range of the parameter A 
depend on the particular scalar tensor theory. These solutions can be generated^ from the 
general relativistic Mars's solution using the solution generating methods developed in 



Ref. m 



Below we consider the explicit form of the general solution for some particular scalar- 
tensor theories. 



Some of the solutions were first obtained by solving the corresponding system of partial differ- 
ential equations for nondiagonal G2 cosmologies 



3 



A. Barker's theory 



Barker's theory is described by the functions = $ and Z($) = (4-3$)/2<l>(<l> - 1). 
In the case of Barker's theory the exphcit forms of the functions and /(A) are: 



<^-\t) = 1 - A cos^ (^a^7 - Itj , (4) 

/(A) = 1 - A (5) 

where the range of A is < A < 1. This range can be extended to < A < 1. For A = 
and A = 1 we obtain the Mars's solution and gravi-dilaton vacuum solution, respectively. 
That is why we consider only < A < 1. It should be noted that the range of the parameter 
A is crucial for the curvature invariants. It is easy to see that the gravi-dilaton vacuum 
solution corresponding to A = 1 has divergent curvature invarinats because of the conformal 
factor = sin^ (aA/7 - li) • 

B. Brans-Dicke theory 

Brans-Dicke theory is described by the functions = $ and Z{^) = uj/^ where lj 

is a constant parameter. Here we consider the case uj > —3/2. The explicit form of the 
functions and /(A) in the Brans-Dicke case is the following: 

$-1/2 (t) = ^exp 

+ (l-A)exp(^-.y^tj. (6) 

/(A)=4A(1-A), (7) 

Here the range of the parameter A is < A < 1. The solution exists for A = and 
A = 1, too. In these cases, however, we obtain a gravi-dilaton vacuum solution which is just 




the Pimentel's solution ||30[. That is the reason we do not consider these limiting values of 
A. The solution is invariant under the trasformations A < — > 1 — A and t < — > —t. In this 
generalized sense, we can say that the solution is even in time. 

C. Theory with "conformal" coupling 

The theory with "conformal" coupling is described by the functions F[^) = 1 — and 
Z($) = 1. In this case we have: 



= 1 + 4A(1 - A) sinh^ \^a^l-ltj , (8) 

/(A) = (l-2Af. (9) 

The range of the parameter AisO<A<l/2. For A = 1/2 we obtain a gravi-dilaton 
vacuum solution which is well-behaved and can be included as a limiting case. 
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D. Z($) = (172 _ 3$)/2$2 theory 



Here we consider the scalar-tensor theory described by the functions = $ and 

Z($) = _ 3$)/2$2 ^i^g^g n > 0. The exphcit forms of and /(A) are: 



= (l + ia^^t) + A, (10) 
/(A) = A. (11) 

Here, the range of the parameter is < A < oo. 

E. Z($) = i($2 - 3$ + 3)/$($ - 1) theory 

The theory with F($) = $ and Z($) = - 3^ + 3)/^{^ - 1) possesses the following 
solution: 

A^ 

^ '^^^ ^ A2 + (l-A2)sin2(Aav/^t)' ^^^^ 
/(A) = A^. (13) 

In order for the dilaton field in this solution to have positive energy we should restrict 
the range of the parameter A to < A < 1. 



Using the solution generating methods developed in Ref. we can generate nondiagonal 
G2 inhomogeneous cosmological solutions in many other scalar-tensor theories different from 
those considered above. However, the solutions we have presented here are expressed in a 
closed analytic form and they are also representative and cover a wide range of the possible 
behaviors of the scalar-tensor solutions which can be generated from Mars's solution. 

Let us consider the main properties of the found solutions. The metric functions, the 
gravitational scalar (the dilaton) and the fluid energy density are everywhere regular. The 
space-times described by our solutions have no big-bang nor any other curvature singularity - 
the curvature invariants Ji = C^i/a/jC^'^"^, I2 = R^uR^^ , and /a = R^ are regular everywhere. 
The solution possesses a two dimensional abelian group of isometries inherited from the seed 
Mars's solution and generated by the Killing vectors d/dz and d/dcf). In addition, the metrics 



have a well defined axis of symmetry and the elementary flatness condition ||3^ is satisfied. 



Since the presented solutions are conformally related to the Mars's solution, the spacetimes 
descibed by them are globally hyperbolic. In fact, the global hyperbolicity can be proved 
independently as a consequence of the poof of the geodesic completeness presented below. 

The existence of two Killing vectors gives rise to two constants of motion along the 
geodesies: 



X 



, N 9*^0 0''^'^ ,dz 9 
cosh(2at)r^-- H (— + ar^ 



ds cosh(2at) ds ds 

F-^{^{t)) dz ^^^2d<P- 
cosh(2at) ds ds' 



(14) 
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The affinely parameterized causal geodesies satisfy 



F-($(t)){e-^^^^osh(2at)[(|)^-(|)^] 



{K - Lar 



2\2 



cosh(2at) 

' " r2F-2($(t))cosh(2at) 



} 



(15) 



where e = and 1 for null and timelike geodesies, respectively. Taking into account ([Hj) 
and (|T5|) the geodesic equations for t and r can be written in the following form: 



d 
ds 



F-i($(t))e^'^'"'cosh(2at) 



dt^ 
ds 



F($(t))e-^'^'"' cosh-^(2at)M(t, r)dtM{t, r) 



d_ 

ds 



F-\(!>{t))e^^ " cosh(2at) 



dr" 
ds 



(16) 



(17) 



where 



-F($(t))e-^"'"' cosh"\2at)M(t, r)drM{t, r), 



M(t,r) = F-^/2($(^))g(l/2)7a2r.2^Qgj^l/2^2at) 



cosh(2at) 



(ir - Lar 



2\2 



nl/2 



F-i(<l>(t)) r2F-i(<l>(t))cosh(2at) 



(18) 



To demonstrate the geodesic completeness of our metric, we have to show that all non- 
spacelike (i.e. causal) geodesies can be extended to arbitrary values of the affine parameter. 
We shall consider only future directed geodesies. The past directed geodesies can be treated 
analogously. 

and 



First we consider null geodesies with K = L = 0. For them we have ^ 



I — I 

I ds I 



4- ( F-\^t))e^'''''\osh(2at)^ 
ds \ ds 



0. 



(19) 



After integrating we obtain 

dt 
ds 



C [F-\^t))e^''''" cosh{2at)) ' 



(20) 



where C > is a constant. Taking into account that for each of our solutions there exists 
a constant B such that 



0<B< F-\<^{t)) 

for arbitrary values of t and fixed parameter A, we obtain 

dt dr C 
ds ds ~ B 



(21) 



(22) 



Therefore the geodesies under consideration are complete. 
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Now let us turn to the general case when at least one of the constants e, K or L is different 
from zero. Here we shall use a method similar to that for diagonal metrics described in Ref. 
p!6| . Let us parameterize ^ and ^ by writing : 



dt _ F($(t))e-^" " 
ds 
dr 
ds 



cosh(2at) 
F(<l>(t))e-^"' 



M(t, r) cosh(i;), 
-M{t, r) sinh(t;). 



cosh(2at) 

Substituting these expressions in the equations for t and r we obtain 

dv F(<l>(t))e-T"^'''' 



ds 



cosh(2at) 



— [dtM{t,r) smh{v) 
+drM{t,r) cosh(i;)] 



or equivalently 



dv 
ds 



2M{t,r 



-{r+(t,r)e" + r_(t,r)e-"} 



where 



r+(t,r) 



, ldln[F-i($(t))] 2 
atanh(2at) + ^ + 7a V 



{K - Lar 



2\2 



r2F-i($(t))cosh(2at) 
cosh(2at) 



dt 



'ya r 

r 



2Lar 



+ 



K — Lar"^ 
2atanh(2at) + ■ya^r 



r_(t,r) 



, , , 1 d\n\F-\^(t))] 2 ■ 
-atanh(2at) - ^ ^ + 7a V 



{K - Lar 



2\2 



r2F-i($(t))cosh(2at) 
cosh(2at) 



2Lar 



'ya r 

r 



K — Lar"^ 
2atanh(2at) + 7a^r 



(23) 
(24) 



(25) 



(26) 



(27) 



In order for the geodesies to be complete 4^ and ^ have to remain finite for finite values 



dr 



of the affine parameter. In fact, it is sufficient to consider only ^, since ^ and ^ are related 
via (|15|). The derivatives ^ and ^ are regular functions of t and r, and the only problem 
we could have appear when r approaches r = for if 7^ 0. We shall show, however, that r 
cannot become zero for if 7^ 0. 

First we consider geodesies with increasing r (i.e. v > 0). In this case it is not difficult 
to see that the term 



Fm))e 



cosh(2at) 



-M{t,r) 



(28) 
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in eqn. (^) can not become singular (for increasing r). Therefore, ^ could become 
singular only for v. We shall show, however, that v can not become singular for finite values 
of the affine parameter. For increasing r, v cannot diverge since for large t (large r ) the 
derivative ^ becomes negative. Indeed, for all exact solutions presented here, there exists 
a constant i?i > such that 

"'■'l^-JfW)' |< (29) 

for arbitrary t and fixed A. 

Therefore, as can be seen from eqns. (^), the terms associated with the constant £, K 
and L are all positive for large values of t.As a consequence we obtain that the functions 
r+(t, r) and r_(t, r) are positive, i.e. ^ < for large t(large r).0 

In the second case,when r decreases (t; < 0), the problem comes from r = when K ^ Q. 
The geodesies with = can reach the axis r = without problems and then continue with 
r > (t; > 0). When K ^ Q^v cannot diverge for finite values of the affine parameter. This 



dr 
ds 



follows from the fact that the derivative becomes positive for small r (large t) as can be 
seen from Eqns. (|27|) and (^6]), taking into account that r_|_(t, r) is exponentially suppressed 
compared with r_(t,r). The positiveness of the derivative ^ when the geodesies are close 
to the axis r = prevents the radial coordinate from collapsing too quickly and reaching the 
axis. The fact that r can not become zero for K ^ Q may be seen more explicitly as follows. 
When r approaches zero the dominant term is that associated with K and the other terms 
can be ignored. So, for small r the geodesies behaves as null geodesies with L = 0: 

*-^(*('»'"^°"''M(.)co=h(.), (30) 



ds cosh(2at) 
dr F($(t))e-^' 



ds cosh(2at) 
dv F(<l>(t))e-T«''- 



-M(r) sinh(i;), (31) 
drM{r) coshiv) (32) 



ds cosh(2at) 
where M(r) = lMg(i/2)7a^r-^_ Hence, we obtain the orbit equation 

dr Mir) , , , 

_ = - \/ tanh i; . (33) 

Integrating, we have 

g-(l/2)7a2r-2^ = Ci COsh(i;) (34) 

where Ci > is a constant. Since cosh(i;) > 1, r can not become zero. 
From the proof of the geodesic completeness it follows that every maximally extended null 
geodesic intersects any of the hypersurfaces t = const. According to this a sufficient 



^In fact, the function r_(t,r) is exponentially small compared with r+(t,r) and may not be 
considered. 
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condition that the hypersurfaces t = const are global Cauchy surfaces. Therefore, the 
solutions are globally hyperbolic. 

We have explicitly proven the geodesic completeness of the solutions using their particular 
properties. The geodesic completeness can be proved independently by considering the solu- 
tions from a more general point of view. In Ref. |^(see also Ref. [^), Fernandez- Jambrina 
presented a general theorem providing wide sufficient conditions for an orthogonally transi- 
tive cylindrical space-time to be geodesically complete. It can be verified that the solutions 
presented here satisfy all conditions in the Fernandez- Jambrina's theorem and therefore they 
are geodesically complete. 

New diagonal solutions can be obtained from (0) as a limiting case. Taking a ^ 
and keeping a^7 = (3 fixed, we obtain the following diagonal inhomogeneous cosmological 
scalar-tensor solutions: 

ds^ = F-\^t)) [e^'^'i-dt^ + dr^) + r^d(P' + dz^ 

8nG,p = Pf{\)e-^''"F%^t)), (35) 
= F-i/2($(t))e(i/2)/3r2 5j 

where ay/j — 1 should be replaced by ^/J3 in the explicit formulas for 
We have proven that the solutions presented in the present paper are geodesically com- 
plete. This result is not in contradiction with the well-known singularity theorems because 
in our case the strong energy condition is violated in the Jordan frame. This can be explicitly 
seen by calculating the components of the Ricci tensor. All components are bounded except 
for Rtr = —r'ja'^dtln{F[^{t)]}. Therefore, for large enough r, one can always find timelike 
and null vectors v'^ such that R^^v^v'^ < i.e. the strong energy condition is violated. How- 
ever, the situation is different in the Einstein frame. The Einstein frame metric gj^^ is just 
the Mars's metric and it is geodesically complete as we have already mentioned. Since the 
energy conditions are satisfied in the Einstein frame it remains to see which other conditions 
of the singularity theorems are violated. The space-time described by the metric gj^^, does 
not contain closed trapped surfaces. In order to prove this we will employ the techniques 



of differential geometry described in Refs. |]39[ and [§0|- Let us consider a closed spacelike 



surface S and suppose that it is trapped. Since the surface is compact it must have a point 
q where r reaches its maximum. Let us denote Vmax = R on a constant time hypersurface 
t = T. For the traces of both null second fundamental forms at q, it can be shown that (see 



Refs. m and |0|) 

g-(l/2)7a2K2 
K^\g > . . , _. > 



V2Rcosh^/\2aT) 



-(l/2)7a2ij2 



Ks\q < F < 0. (36) 

The traces have opposite signs so that there are no trapped surfaces. 

Our solutions are stiff perfect fluid cosmologies and, therefore, the natural question which 
arises is what happens if the fluid is not stiff. In this case, however, the situation is much 
more complicated. In contrary to the stiff fluid case, the dilaton-matter sector does not 
posses nontrivial symmetries which allow us to generate new solutions from known ones. The 
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only way to find exact solutions is to attack directly the corresponding system of coupled 
partial differential equations. This question is currently under investigation. 

Summarizing, in this work we have presented new nondiagonal G2 inhomogeneous stiff 
perfect fluid cosmological solutions in a wide range of scalar-tensor theories. The found 
solutions have no big-bang nor any other curvature singularity. The gravitational scalar 
(dilaton) and fluid energy density (pressure) are regular everywhere, too. Moreover, the 
solutions are globally hyperbolic and geodesically complete. To the best of our knowledge, 
these solutions are the first examples of nonsingular G2 inhomogeneous perfect fiuid scalar- 
tensor cosmologies with a nondiagonal metric. 

I would like to thank V. Rizov for discussions and especially L. Fernandez- Jambrina for 
his valuable comments on the geodesic completeness of the orthogonally transitive cylindrical 
spacetimes. My thanks also go to J. Senovilla for sending me some valuable papers. This 
work was supported in part by Sofia University Grant No 459/2001. 
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